Significant X-ray Line Emission in the 5-6 keV band of NGC 4051 
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ABSTRACT 

A Suzaku X-ray observation of NGC 4051 taken during 2005 Nov reveals line emission at 
5.44 keV in the rest- frame of the galaxy which does not have an obvious origin in known rest- 
frame atomic transitions. The improvement to the fit statistic when this line is accounted for 
establishes its reality at > 99.9% confidence: we have also verified that the line is detected in the 
three XIS units independently. Comparison between the data and Monte Carlo simulations shows 
that the probability of the line being a statistical fiuctuation is p < 3.3 x 10~^. Consideration 
of three independent line detections in Suzaku data taken at different epochs yields a probability 
p < 3 X 10""'^^ and thus conclusively demonstrates that it cannot be a statistical fiuctuation in 
the data. The new line and a strong component of Fe Ka emission from neutral material are 
prominent when the source fiux is low, during 2005. Spectra from 2008 show evidence for a line 
consistent with having the same flux and energy as that observed during 2005, but inconsistent 
with having a constant equivalent width against the observed continuum. The stability of the 
line flux and energy suggests that it may not arise in transient hotspots, as has been suggested for 
similar lines in other sources, but could arise from a special location in the reprocessor, such as 
the inner edge of the accretion disk. Alternatively, the line energy may be explained by spallation 
of Fe into Cr, as discussed in a companion paper. 
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1. Introduction 

New X-ray data from high-throughput missions 
such as XMM-Newton and Suzaku, along with the 
high spectral resolution afforded by the Chandra 
High Energy Transmission Grating (HETG) has 
led to the discovery of several important spectral 
signatures in Active Galactic Nuclei (AGN) that 
were not detectable using previous missions. 

A joint XM M-Newton/ Chandr a observation 
of NGC 3516 (|Turner et al. 2002[ ) revealed the 
first detection of narrow line emission at 5.6 and 
6.2 keV (the latter was able to be resolved from 
the strong Fe Ka line owing to the relatively high 
spectral resolution afforded by the HEG grat- 
ing). While the line energies matched those ex- 
pected for ionized species of Cr and Mn, the line 
strengths exceeded what might be expected from 
illumination of material with cosmic abundance 
ratios. One explanation considered was enhance- 
ment of Cr and Mn abundances from spallation 
of Fe. However, this was initially disfavored as 
the line r atios did not agree well with those pre- 
dicted by ISkibo ('19971) for spallation of disk gas 
([Turner et al. 20021) . The alternative possibility, 
of Doppler-shifted components of Fe emission, 
led to a suggestion of emitting hotspots on the 
accretion disk surface - with an expectation of 
observable shifts in line flux and energy as the 
hotspot traverses its orbit. Further examples 
of the phenomenon, dubbed 'transient Fe lines', 
soon came from observations of other AGN, with 
many lines reported in the 5-6 keV r egime (e.g. 
Yaqoob et al. 20031 : iTurner et al. 20ol . 

In this paper we report the highly significant 
detection of an emission line at 5.44 keV in the 
rest-frame of the nearby narrow line Seyfcrt 1 
type AGN NGC 4051. The systemic redshift of 
the host galaxy is 2; = 0.0023 that yields a dis- 
tance 9.3Mpc for Hq — 74kms~^Mpc~^ assum- 
ing the redshift to arise from the Hubble flow. 
However, for such nearby galaxies a more reli- 
able estimation of distance can be obtained from 
use of the T uUy-Fisher rel ation, giving a distance 
of 15.2 Mpc ( Russell 2004f l in this case, which we 
adopt in this paper. 

NGC 4051 was observed by Suzaku in 2005 and 
2008. The 2005 observati ons have previously been 
described and analyzed by lTerashima et aL (;2009h . 
who showed that the X-ray source was highly vari- 



able on short timescales, and that a spectral model 
including variable partial-covering absorption was 
required to fit the data. Here we make a joint 
analysis of the 2005 observation together with the 
new 2008 data in which we concentrate specifi- 
cally on the detection of the narrow emission line 
and its variability. We then discuss the possible 
origin of the line in the context of several popular 
models. 

2. Observations 

The Suzaku X-ra y Imaging Spectrometer (XIS 
Kovama et al. 2007h instrument com prises four X- 
ray telescopes (jMitsuda et al. 20071 ) each with a 
CCD in the focal plane. XIS CCDs 0,2,3 are con- 
figured to be front-illuminated (FI) and provide 
useful data over ~ 0.6 — 10.0 kcV with energy res- 
olution FWHM - 150 eV at 6keV. XIS 1 is a 
back-illuminated CCD and has an enhanced soft- 
band response (down to 0.2 keV) but lower area at 
6 keV than the FI CCDs as well as a larger back- 
ground level at high energies, consequently this 
detector was not used in our analysis. Suzaku also 
carries a non-imagin g, collimat ed Hard X-ray De- 
tector (HXD Takahashi et al. 2007 ) whose PIN in- 
strument provides useful data over 15-70 keV for 
bright AGN. 

Our analysis used all available Suzaku observa- 
tions of NGC 4051 comprising data from 2005 Nov 
10-13 (observation identifier 700004010) and 2008 
Nov 6-12 (703023010) and 23-25 (703023020) as 
summarized in Table 1 . The data were reduced us- 
ing v6.4.1 of HEAsOFT and screened to exclude: 
i) periods during and within 500 seconds of the 
South Atlantic Anomaly (SAA) , ii) with an Earth 
elevation angle less than 10° and iii) with cut-off 
rigidity > 6 GeV. The source was observed at the 
nominal center position for the XIS during 2008 
and at the nominal center position for the HXD 
during 2005. The FI CCDs were in 3 x 3 and 5 x 5 
edit-modes, with normal clocking mode. For the 
XIS we selected events with grades 0,2,3,4, and 
6 and removed hot and flickering pixels using the 
SISCLEAN script. The spaced-row charge injec- 
tion (SCI) was used. The XIS products were ex- 
tracted from circular regions of 2.9' radius with 
background spectra from a region of the same size, 
offset from the source (avoiding the calibration 
sources at the edges of the chips). The response 
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and ancillary response files were created using xis- 
RMFGEN v2007 May and xissimarfgen v2008 
Mar. 

NGC 4051 is too faint to be detected in the 
HXD GSO instrument, but was detectable in the 
PIN. For th e analysis we used t he model "D" 
background ( Fukazawa et al. 2009() . As the PIN 
background rate is strongly variable around the 
orbit, we first selected source data to discard 
events within 500s of an SAA passage, we also 
rejected events with day/night elevation angles 
> 5°. The time filter resulting from the screen- 
ing was then applied to the background events 
model file to give PIN model-background data 
for the same time intervals covered by the on- 
source data. As the background events file was 
generated using ten times the actual background 
count rate, an adjustment to the background 
spectrum was applied to account for this factor. 
HXDDTCOR v2007 MAY was run to apply the 
deadtime correction to the source spectrum. To 
take into account th e cosmic X-ray bac kground 
(|Boldt fc Leiter 1987t iGruber et al. 19991 ) xspec 
v 11.3.2ag was used to generate a spectrum from 
a CXB model ( Gruber et al. 19991) normalized to 
the 34' X 34' Suzaku PIN field of view, and com- 
bined with the PIN instrument background file to 
create a total background file. The mean exposure 
times for XIS are given in Table 1. The exposure 
times for the PIN were 112, 204 and 59 ks for 2005 
Nov 10, 2008 Nov 6 and 2008 Nov 23, respectively. 

Spectral fits used data from XIS 0, 2 and 3, 
for the 2005 data. As use of XIS2 was discon- 
tinued after a charge leak was discovered in Nov 
2006, the 2008 analysis used only XIS and 3. In 
this paper, XIS data were fit over 3.0 — lOkeV. 
PIN data were fit simultaneously with the XIS, in 
the range 15-50 keV. In the spectral analysis, the 
PIN model fiux was increased by a factor 1.16 for 
2008 data and 1.18 for 2005 data, which are the 
appropriate adjustments for the instrument cross- 
calibration at those epochs of the observation. XIS 
data were binned at the HWHM instrumental res- 
olution while PIN data were binned to be a min- 
imum of 5(7 above the background level for the 
spectral fitting. 



3. Spectral Fitting 

The mean Suzaku count rate over 0.5-10 keV 
during the low state of 2005 was 0.45 XIS 
counts"^ per FI XIS and 0.04 PIN count s'^ 
During the high-state of 2008 Nov 6 it was 
2.06 XIS count s-i per FI XIS and 0.07 PIN 
count s"'^ (Table 1). The count rates recorded 
correspond to an observed 2-10 keV flux range 
0.87 - 2.42 X 10-"ergcm-2s-i and a 10-50 keV 
fiux range 2.44 - 3.63 x IQ-^ergcm-^s'^ The 
2005 dat a represent a historica lly low state, as 
noted bv lTerashima et al. f2009h and during 2008 
the source was close to the historical average fiux 
level. The 2008 Nov 6 Suzaku observation was 
accompanied by a contemporaneous HETG expo- 
sure that reveals a wealth of emission and absorp- 
tion features from several zones of ionized gas. 
The HETG an alysis will be des cribed in detail in 
another paper (|Lobban et al. 2010.V 



3.1. Fitting the PC A Offset Component 

To assist in modeling the full-band Suzaku data, 
we first performed a decomposition of the events 
from all three Suzaku observations using Princi- 
pal Components Analysis using the SVD method 
and code of Miller et al. ^20071 ). A full descrip- 
tion of the PC A decompositio n of these data is re- 
ported bv lMiller et al. C2009I) . In summary, PC A 
is a mathematical decomposition of the data into 
orthogonal 'eigenvectors'. The dominant varia- 
tions are ascribed to eigenvectors of the lowest 
orders. In the case of NGC 4051 the data were 
found to be well-described by a steady "offset" 
co mponent having a hard spectrum (Figure [T]and 
see iMiller et al. 200 9) while the first order vari- 



able component, eigenvector 1, is consistent with 
an absorbed power-law of constant slope whose 
variations in intensity dominate the spectral vari- 
ability of_yie_souree_(consistent with the analy- 
sis of iTerashima et al. 20091) . The physical ori- 
gin of the offset component is of great interest. 
The low-state spectrum of NGC 4051 is com- 
posed of this component plus some contribution 
by eigenvector 1. Fitting the lowest fiux sub- 
set of the 2005 data, Terashima et al. 20091 found 
the hard component to be best described using a 
partially-covered refiection component. Our anal- 
ysis of the PCA offset component provides an- 
other way to probe the lowest fiux levels of the 
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source and, while we find an acceptable fit using a 
combination of reflection (modeled using PEXRAV) 
plus a contribution from the powerlaw contin- 
uum (r = 2.3, fixed from fitting eigenvector!), 
both emission components require absorption by a 
complex of ionized gas. Fortunately the multiple 
layers of ionized gas can be well constrained us 



ing HETG and LETG data (ICoUinge et al. 2001 



ISteenbrugge et al. 2009t Lobban et al. 2010 ). re- 
ducing the degeneracy of the Suzaku fit. In ad- 
dition to these components, FeKa line emission 
is evident in the offset spectrum at a fitted en- 
ergy E = 6.398 ± 0.023 keV having a width a = 
0.045t^;^^^keV and fiux n = 1.90 ± 0.26 x lO'^ 
photons cm~^ s~^. (Throughout this paper, errors 
are quoted at 90% confidence for the appropriate 
number of interesting parameters in the fit.) The 
resulting fit statistic was — 130/153 d.o./. 

Intriguingly, the normalization of the reflection 
continuum is extremely high relative to that of the 
primary continuum. The inconsistency of reflec- 
tion and continuum strengths may indicate that 
the partial-covering absorption comprises a larger 
part of this component than currently modeled. 
The relatively hig h flux of the hard com ponent 
was also noted by iTerashima et al. r2009h based 
up on fltting the 2005 d ata, and is fully explored 
bv lLobban et al. (2010l ). 

In addition to the prominent narrow Fe K line, 
a strong excess of line emission is evident in the 
PGA oflFset spectrum, at 5.44 keV (Figure [1]). No 
significa nt residua ls appear in eigenvector 1 at that 
energy (jMiller et al. 2009 '). The new line is thus 
consistent with having an origin associated with 
the hard offset component and neutral component 
of Fe Ka emission. 

3.2. Fitting the 2008 Contemporaneous 
Chandra and Suzaku data 

As the 2008 Nov 6 epoch provides contempora- 
neous HEG data for the Suzaku observation, we fit 
those spectra together over 3-8 keV, combining the 
superior energy-resolution of HEG and the good 
statistical quality of the FI XIS units (0,3) to ob- 
tain the best possible constraints on the width and 
energy of the neutral component of Fe Ka. A joint 
fit to those data using a simple powerlaw plus sin- 
gle Gaussian line showed residual excesses in the 
5-6 keV band in Suzaku XIS data (Figure [2]). 



We reflt, adding to the model a single com- 
ponent of absorption, a Gaussian representa- 
tion of the K/3 line fixed at a rest-energy of 
7.05 keV with a flux li nked to be 13.5% of the 
Kg component ( e.g. Leahv fc Greighton 19931 : 
Palmeri et al. 20031 ). a narrow absorption line 
detected at a n observed energy of T.lkeV (as 
discovered by [Pounds et al. 20041 and confirmed 



by ITerashima et a l. 2009) plus a component to 
model the narro w line emission eviden t at 6.62 keV 
in HEG data (jLobban et al. 2010l ). The flt 
yielded a line energy E = 6.410 ± 0.015 keV, 
width (T = 50^33 eV and normalization n = 
1.90 ± 0.395 X lO-^photonscm-^s-^ The equiv- 
alent width for the Fe Ka line component is 95 eV 
against the total continuum level observed during 
2008. The constraint on the width of the Fe Ka 
line is equivalent to a FWHM velocity-broadening 
of 5540^3554 km s^^ (90 percent confidence inter- 
val) . Assuming the true line centroid energy to be 
6.40 keV, the 90% confldence constraint on energy 
from the combined HEG/XIS data limits the bulk 
velocity of this emitter to 250 ^ v > — 1200 km s"""^ 
where negative velocity denotes outflow. 

As the PGA decomposition indicates the pres- 
ence of a we ak broad compone nt of Fe Ka on 
eigenvector 1 ( Miller et al. 2009[ ). we tried adding 
a second component of Fe Ka emission to the 
model, with the energies of both line compo- 
nents linked. This addition improved the flt- 
statistic by A^^ — 5, giving the same line en- 
ergy as for the single component model, with 
(Ti = 1^^'^ eV and normalization ni ~ 1.12 ± 
0.32 X 10-5photonscm-2s-i (EW= 41 eV), era = 

= 1.07 ±0.45 X 
eV). Although 



145_52 and normalization n2 ■ 
lO-^photonscm-^s-i (EW= 43 
the two-component parameterization of Fe Ka of- 
fers only a marginal improvement to the flt, it is of 
interest with regard to consistency with the PCA 
decomposition, which indicates the weak broad 
component of Fe Ka emission to be present on 
eigenvector 1. 

In subsequent fits we flx the Fe Ka model line 
width at (T = 50 eV obtained from the simple flt, 
while bearing in mind the possible more complex 
solution for the line proflle. As the PCA is consis- 
tent with a common origin for Fe Ka emission and 
the line at 5.44 keV, and as the widths of the two 
lines are consistent (and as there is no other use- 
ful constraint available for the indeterminate line 
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width) hereafter we fixed both to a ~ 50 eV in 
spectral fitting. 

3.3. Fitting the 2005 Suzaku data 

To assess the significance of the line at 5.44 
keV we returned to direct fitting of the 2005 data, 
where the source was observed to be at a low flux 
and where PCA indicates that the line would have 
the highest equivalent width. We fit the summed 
2005 data from the FI XIS units 0, 2 and 3 using 
an absorbed continuum model, applying it to the 
3-10 keV band for the purpose of examining the 
data for features of interest. A Gaussian line was 
included in the fit, to account for the Fe Ka emis- 
sion, fixed at the energy and width found using 
HEG/XIS. As evident in Figure [3] the mean spec- 
trum from 2005 data shows an excess of counts at 
5.44 keV with additional structure evident around 
6 keV. We performed a more complex fit to prop- 
erly assess the line strengths and significance. We 
added to the model a Gaussian representation of 
the Fe K(3 line fixed at a rest-energy of 7.05 keV 
with a flux linked to be 13.5% of the Ka com- 
ponent. The model also included the narrow ab- 
sorption line detected at an observed energ y of 
7.1keV as di scovered bvlPounds et al. f2004l) and 
confirmed bv lTerashima et al. (2009 ). plus a com- 
ponent to m odel the narrow line emission evident 
at 6.62 keV (|Lobban et al. 2010l ). The fit yielded 
= 176/101 dof and the strongest unmodeled 
feature remains evident as an excess of counts at 
5.44 keV. 

Addition of a Gaussian line (cr — 50 eV) re- 
duced the fit-statistic by Ax^ — 32, yielding 



E = 5.44 ± 0.03 keV, flux 
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photons cm~ s~ and equivalent width 46 ± 16 
eV. Addition of a second line of the same width 
yielded a further improvement Ax^ = 34, E — 
5.95±0.05keV, line flux n = 5.0511;^^ x 10"*^ pho- 
tons cm~^s~"'^, with an equivalent width of 44^i[g 
eV and a final fit statistic x^ — 117/97 d.o./. In 
this fit the flux and equivalent width of the Fe Ka 
line were n = 1.59 ± 0.23 x 10"^ and 195 ± 24 eV 
respectively. 

To further examine the source behavior at low 
flux levels we isolated the lowest flux subset of the 
2005 data using an intensity fllter based upon the 
0.5-10 keV band count rate, taking data below a 
threshold of 0.47 ct s'^ (per FI XIS). The intensity 



cut effectively removed the periods of source flar- 
ing from consideration. This filtering thus resulted 
in a very-low-state spectrum that had a mean 2- 
lOkeV flux 6.1 X lO^^^crgcm^^s^^. Again, fit- 
ting a simple absorbed continuum model over 3- 
lOkeV leaves two strong positive residual features 
between 5-6 keV (Figure |4]) confirming the promi- 
nence of the 5.44 keV line at low flux levels and 
indicating the possible presence of emission at 5.95 
keV. 

3.4. The Full Model 

To assess the full properties of the source and 
investigate the robustness of the line detections to 
the continuum model used, we then fit over 0.75- 
50 keV (the full band of the Suzaku data). We 
simultaneously fit spectra from the three Suzaku 
observations. The 2008 HETG exposure (that 
overlapped the 2008 Nov 6 Suzaku observation) 
shows a wealth of absorption and emission lines 
that can be modeled using three zones of ion- 



ized gas (jLobban et al. 20101 ): those soft-band ab- 
sorbers were fixed in the Suzaku fit, after which 
two additional gas layers were required to accu- 
rately model the source across the full Suzaku 
bandpass. The fit also included an Fe Ka emis- 
sion line fixed at a width of ct — 50 eV as 
before. Finally, an absorbed ionized refiector 
was included to account for the curvature in the 
soft band. The full model is detailed in Ta- 
ble 2 and the fit is illustrated in Figure [5] The 
marked spectral variability observed (Figure [5]) 
can be accounted for by allowing variations in 
the covering fraction of one absorbing layer (hav- 
ing Nh 1023cm-2,log^ - 0.18), a solution 
that is si milar to other well-studied AGN of this 
class ( e.g. IPounds et al. 200l iRisaliti et al. 2007 



iMiller et al. 2008t iTurner et al. 200*8^ This fit 
yielded x^ = 947/547 d.o.f. with the overall curva- 
ture modeled well across all epochs and the resid- 
ual contributed mainly by some unmodeled 
spectral features. The addition of a line (ct = 50 
eV) constrained to lie in the range 3-8 keV yielded 
a detection at 5.44±0.03 keV and improved the fit 
statistic by Ax^ — 20; addition of another line un- 
der the same constraints gave 5.98±0.05 keV and 
Ax^ — 9 improvement to the fit. We conclude that 
the high significance for the detection of the 5.44 
keV line is robust to the continuum used while the 
significance of the detection of a line at 5.95 keV 
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is sensitive to the continuum form assumed. 

3.5. A Critical Examination of the Reality 
of the Lines 

Before proceeding any further with the mod- 
eling, we performed several checks to be assured 
that the new line is attributable to the AGN and 
is significant. First, to confirm that the line is not 
an artifact of poor background subtraction, we ex- 
amined the background spectrum and calibration 
source data. First we note that the background 
comprises just 2.5% of the total count rate in the 
5-7 keV band for the XIS spectral data when the 
source is at its lowest flux level, during 2005 and 
1.3% when the source is brighter, during 2008. Re- 
garding the line at 5.44 keV, we found there to be 
no line emission evident at that energy in the back- 
ground spectrum. Fitting the background spectra 
from the three Suzaku observations we obtained 
an upper limit (90% confidence) for the flux of 
a line at 5.44 keV in the background spectrum 
n < 4.14 X 10"® photons cm"^s"\ i.e. < 1% of 
the detected line flux. 

The lack of a feature at comparable flux or 
equivalent width in the background data also rules 
out an origin of the 5.44 keV line as a detector fea- 
ture. We then examined each XIS independently, 
and found the 5.44 keV line to be significantly de- 
tected in each XIS unit independently, further ver- 
ifying the reality of the reported feature (Figure 

El). 

Regarding the more tentative line at 5.95 keV: 
the XIS chip does include a Mn calibration source 
emitting a line at 5.9 keV for the purpose of cal- 
ibration of the XIS energy scale. The calibration 
source is located at the chip edge and the counts 
from that source are obviously conservatively ex- 
cluded from any source and background extraction 
cells for scientific analysis. However, we do detect 
a weak Mn Ka emission in the background spec- 
trum. The Mn contamination is at a level < 10% 
of the measured line strength in the AGN spec- 
trum. As the source extraction cell is further from 
the Mn calibration source than the background 
cell and therefore subject to less contamination, 
we estimate the Mn line contamination from the 
calibration source to be < 10% for the 2005 spec- 
trum. In spectral fitting, the greatest concern for 
the Mn line is the possible contamination by the 
weak broadened Fe Ka emission evident in the 



fit to eigenvector 1 ( Miller et al. 20091 ) which may 
lead to the strength of the Mn line being over- 
estimated from simple fits. Because of the various 
issues associated with a clean measurement of any 
line at 5.95 keV and the sensitivity to the con- 
tinuum form assumed we concentrate only on the 
strong line detection at 5.44 keV. 

To confirm the significance of the line we 
performed Monte Carlo simulations us ing the 
method described in iPorquet et al. (20041 ) and in 
Markowitz et al. (2006ll We took the null hy- 



pothesis to be that the spectrum is simply an 
absorbed power-law continuum with parameters 
derived from fitting the broad-band data but al- 
lowing for the statistical uncertainty on the con- 
tinuum parameters and including the narrow Fe 
Ka line whose prese nce is w ell-established in this 
source (Lobban et al. 2010h . We used the xspec 
command fakeit to create 3000 fake Suzaku spec- 
tra with photon statistics expected from the 2005 
exposure, assuming the same instruments to be 
operational as for the actual observation. The sim- 
ulated data were grouped to the HWHM energy- 
resolution of the instruments, the same as the 
observational data. Following the procedure used 
to test the real data for the presence of a narrow 
line, we fitted each fake spectrum to obtain the 
values of Ax^ obtained from statistical fluctua- 
tions in the data. To map the distribution of Ax^ 
across the simulated spectra, each simulated spec- 
trum was fitted over the 3-10 keV energy range, 
stepping through using energy bins whose centers 
were increased in increments of 100 eV. The line 
energy was allowed to be free within each energy 
bin tested and the value of Ax^ was recorded at 
each point in the spectrum. This method makes 
no assumptions about the energy at which a line 
might be detected and over the course of the test- 
ing, all energies are tested for the presence of a 
line. When we fit the simulated data we are testing 
whether we can produce, from statistical fluctua- 
tions, a contribution to ^it the same or greater 
level as found in the actual data at any energy in 
the range of interest. The fits to the simulated 
data yield a distribution of Ax^ for comparison 
with the actual data. We found the most extreme 
statistical fluctuation to yield a Ax^ contribution 
of 19.3 in the 3-10 keV band in the set of simulated 
data (i.e. in 3000 simulations no false line appears 
at any energy contributing Ax^ = 32 as found in 
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2005 data). Thus the probabihty of satisiying the 
null hypothesis is p < 3.3 x lO""* for the line found 
at 5.44 keV. 

3.6. Application of a Disk Hotspot Model 

As an alternative to modeling using individual 
Gaussian lines, we fitted the data using a single 
disk line from a narrow annulus, such that the red 
horn of such a line might explain the peak at 5.44 
keV. We assumed the system to contain a non- 
rotating black hole, that the line is Fe Ka emission 
from neutral material (at 6.4 keV), that the emis- 
sivity pattern across the disk can be described by 
r~'^ where the emissivity index q=-2.5 and that 
the hotspot exists over a narrow annulus of width 
Ar = Ivg. We used the same baseline model as 
for testing the Gaussian lines (section [3^ . 

To fit the 5.44 keV peak as the red Doppler 
horn of a disk-line requires an emitting radius 
21±3rg (where Vg = GMbh/c^ denotes gravi- 
tational radii) in a low inclination system with 
27I3 '^ith — 128/95 d.o./.; the corresponding 
blue horn is then predicted to lie at 6.62 keV and 
the data are consistent with that model. While 
there is evidence for line emission in the 6.5-7 keV 
regime, such lines are commonly observed as emis- 
sion from ionized species of Fe, and thus the iden- 
tification of emission blue- ward of 6.4 keV is cur- 
rently ambiguous. 

3.7. Examination of the Line Variability 

Fits to the mean spectra from each observation 
have shown consistent line fluxes for Fe Ka and 
the line at 5.44 keV and these show an equivalent 
width that appears to have changed as the source 
flux varied. Figure [7] shows the ratio of data in 
the Fe Ka regime to a common local continuum 
fit; the fact that the source spectrum is steeper at 
high flux is reflected in the systematics of the resid- 
uals. To confirm the significance of the change in 
equivalent width we fit the 2008 data with a model 
that fixed the line equivalent width for the feature 
at 5.44 keV to that found during 2005; after refit- 
ting this resulted in a worse fit with Ax^ = 142. 
Alternatively, fixing the fiux of the new lines at 
the 2005 values and refitting, we found Ax^ = 0, 
indicating that the line fluxes may be consistent 
with lines of constant fiux across the data. 

As the limits on line flux provide the potential 



to distinguish between models, we tested all avail- 
able high-quality data in a self-consistent manner. 
Taking the model from section l373l we flxed the line 
energy and width to speciflcally test the constancy 
of the 5.44 keV line. In addition to testing the 
mean spectrum from each Suzaku observation, we 
sub-divided the 2005 and 2008 exposures to sam- 
ple the line more finely in flux. For 2005 an inten- 
sity selection was made on the 0.5-10 keV count 
rate at 0.47 ct s~^ per XIS as before. For 2008 
intensity selections were made < 2.0 ct s~^/XIS 
(ii), 2.0 - 3.0 ct s-i/XIS {i2) and > 3.0 ct s" VXIS 
{i3). The results are shown in Table 3 where both 
the mean fits for each observation are tabulated, 
as well as the intensity-selected results. As noted 
previously, the line is required at a high level of 
confidence in the 2008 Nov 6 data (which is more 
sensitive to the features than the Nov 23 data, ow- 
ing to the long exposure and high flux state of the 
source). 

Table 3 shows line flux along with improve- 
ments in the fit-statistic and equivalent width for 
each fit. Note that the mean observation fits in 
Table 3 show slightly tighter constraints in line 
fiux for 2005 compared to the values noted in sec- 
tion 13.31 because the initial fits to 2005 data had 
the line energy left free. The data are consistent 
with the 5.44 keV line existing at the same fiux 
level throughout the Suzaku observations consid- 
ered in Table 3 and the high significance of the line 
detection across several time slices of data provides 
compelling evidence for the reality of the line; con- 
clusively ruling out the possibility of the line detec- 
tion being a statistical fiuctuation in the spectral 
data. Considering the three independent Suzaku 
detections of the line, as shown in Table 3 Suzaku 
observations 1,2 and 3 yield improvements to the 
fit Ax^ = 32, 19, 19 respectively and a probability 
of all three detections being false isp<3xl0~^^. 

We repeated the fits with the line energy al- 
lowed to be free and found the fitted energy to 
be consistent with 5.44 keV and that the line fiux 
had not been significantly affected by freezing the 
energy. 

To extend the test for line flux variations we re- 
duced and fit the XMM-Newton spectra obtained 
during 2001 and 2002. We followed the standard 
reduction method for the pn data, as detailed by 
Ponti et al. (20061 ) and found the data to be con- 
sistent with the presence of a line at the same fiux 
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as found using Su zaku. An independent a nalysis of 
the XMMAaXa bv lde Marco et al. (20091 ) reported 
the presence of an excess of counts in the 5.4-6.2 
keV band for NGC 4051, compared to their param- 
eterization of the local continuum. Examination 
of the XMM data shows an excess of emission at 
~ 6 keV in the XMM spectra and so tentatively 
supports the possibility of that additional energy- 
shifted line in this source. 

BeppoSax also observed NGC 4051, finding the 
source to be in a very low flux stat e during 1998 
as reported by ICuainazzi et al. (1998,) . We ex- 
tracted and fit the archived BeppoSax Medium En- 
ergy Concentrator spectrum from units 2 and 3 
(combined) and tested for the presence of the line. 
However, the data did not yield a significant de- 
tection of and the upper limits on line flux were 
very loose ~ 4 x 10~^photons cm~^s~^ for a line 
at 5.44 keV. Given the poor constraint obtained 
we do not consider those data any further. 

Combining the results from Suzaku and XMM 
observations (specifically, Table 3 hues 1,2,4,5,7-9 
& 10) and comparing the data to a constant model 
yields — 4.7/7 d.o./. Further to this test, we 
split the data by time instead of intensity and re- 
peated the test, reaching the same conclusion, i.e. 
that the data are consistent with a line of constant 
flux over a timescale of several years and over large 
changes in observed continuum flux. The limits 
on measured line fluxes mean we can rule out line 
variability greater than a factor ^ 2 in the line 
flux sampled on these timescales, across the base- 
line time period considered. 

4. Discussion 

While numerous claims exist in the literature 
for emission lines at unexpected energies (see 
Turner fc Miller 2009aL for a review), the real- 
ity of the lines has bee n questioned by some. 
Vaughan fc Uttlev f 20081 ) considered 38 published 
results on transient emission and absorption lines 
reported in the literature: those authors find 
a linear relationship betwe en the fitted feature 
strength and its uncertainty. Vaughan fc Uttlev (2008[ ) 
noted that observations with more signal appar- 
ently reveal weak lines but do not show tightly 
constrained strong lines which should sometimes 
also show up by chance. The conclusion of 
the I Vaughan fc Uttlev (20081 ) literature review 



was that there is a publication bias in report- 
ing of these results, and that many of the re- 
ported detections are merely statistical fluctu- 
ations. This question has now been addressed 
with two systematic analyses of samples of AGN. 
Tombesi et al. (^2010[ ) present a study of the oc- 



currence and reality of energy-shifted absorp- 
tion lines in a sample of AGN finding a devia- 
tion from the linear relationship of line equiv- 
alent width (EW) and uncertainty found by 



Vaughan fc Uttlev (20081 ) in the sense that their 



distribution showed more significant detections 
of lines in sources studied. Tombesi et al. (2010( ) 
compare their absorption line measurements with 
those of the Fe Ka emission in the same sources 
and show that these two sets of measurements fol- 
low the same distribution in the EW/ uncertainty 
plane. iTombesi et al. i2oi3) conclude that the 
absorption line detections are generally not sta- 
tistical fiuctuations and that the absence of well- 
constrained detections o f strong features in the 



Vaughan fc Uttlev (20081 ) analysis may be due in 



part to a limit on the ability of current X-ray 
instruments to detect such lines. It also would ap- 
pear likely that long observations of bright sources 
may not have been proposed or approved early in 
the XMM and Suzaku missions and so the sample 
of observations completed to date is likely biased 
against those that would have shown tightly con- 
strained energy-shifted lines: such a bias may 
explain at least some of the effect discussed by 
I Vaughan fc Uttlev C2008 ). Another recent study 
bv lde Marco et al. ('20091 ) undertook a systematic 
analysis of a sample of bright Seyfert 1 galaxies 
and confirmed many of the individual detections 
of energy-shifted emission lines claimed in the lit- 
erature, supporting the general reality of the phe- 
nomenon by consideration of the statistics of the 
sample results as a whole. With conflicting views 
in the literature it is clear that the detection of 
significant new examples of the energy-shifted line 
phenomenon is very important at this time. 

4.1. Hard X-ray Line Emission in NGC 
4051 

NGC 4051 can be modeled using a powerlaw 
continuum covered by multiple zones of gas, sev- 
eral in the column density range 10^^ — lO^^'cm"^, 
that impart emission and absorption features to 
the X-ray spectrum. The marked spectral vari- 
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ability with observed flux can be explained by 
changes in covering of the powerlaw continuum 
by one of the high-column a,bsorbers, as fo und 
for other similar AGN (e.g. Pounds et al. 2004 



minating continuum, can be used to set limits on 
the reprocessing gas in which it arises. Following 
Yaqoob et al. and correcting for the con- 



iRisaliti et al. 20071: iMiller et al. 20081: iTurner et al. 
Significant line emission has been observed at 
5.44 kcV in the 2005 Suzaku observation of NGC 
4051. Comparison of low and high-state X- 
ray data for NGC 4051 shows that the newly- 
discovered line appears prominent in the low-flux 
state along with the narrow component of Fe Ka 
emission from neutral g expected if the ob- 

served low state is simply those times when a 
relatively large fraction of the continuum is sup- 
pressed by absorption. 

The most recent measurements indicate NGC 
4051 to have a black hole mass M bh = l-73lo 52 x 
1O^M0 and a radius for the H/3 broad line 
region (BLR) Rr,r p = 1.87t°-^^ light days 
(jPennev et al. 20091). The FWHM in the rms 
spectrum of iDennev et al.l is 1034±41kms^^ (al- 



though those authors used the velocity dispersion 
rather than the FWHM in their mass estimate) . If 
we assume the same geometrical correction factor 
between line width and circular velocity as those 
authors to scale the respective FWHM measure- 
ments, the line width of Fe Ka is indicative of 
an origin at r ~ 1.87(1034/5540)2 ~ 0.065 light 
days, or 2 X 10^^ cm, with a 90 percent confidence 
range 8.6 x 10"^^ — 1.7 x 10^^ cm. As we have 
scaled to the optical reverberation results, this 
radius estimate is secure provided the X-ray and 
optical line-emitting regions have similar struc- 
ture and orientation with respect to the observer, 
but the large uncertainty is dominated by the un- 
certainty in the line width measurement. Further 
to the measurement error is the uncertainty as to 
whether, as suggested by the form of eigenvector 
1, the Fe Ka emission has both broad and nar- 
row components, in which case we would have to 
conclude that we are seeing contributions from re- 
gions both within and outside of the region noted 
above. Given the limited signal-to-noise in the 
regime of Fe Ka in the HEG data this question 
remains open with current data. 

In the single component model for Fe Ka emis- 
sion, spectral fitting to the HEG plus Suzaku 
data yields an Fe Ka emission line of flux n = 
1.90±0.40x 10^5 photons cm^^ g-i. The strength 
of the Fe Ka line emission, normalized to the illu- 



tinuum slope found here, we estimate an efficiency 
2008f ) . for the production of line photons (defined as the 
ratio of line flux to incident flux above the ioniza- 
tion edge) to be x^^kc ^ 0.018. The line measured 
here sets a lower limit on the column density of 
the emitting region Nh > lO^'^cm"^ and in the 
toroidal reprocessor model suggests a global cov- 
ering factor ~ 0.9 with an approximately face-on 
view down the pole of the structure. As it is diffi- 
cult to determine the intrinsic continuum strength 
from the observed continuum the interpretation of 
these data in the context of the toroidal model is 
subject to some uncertainty that in turn, leaves 
the derived global covering factor uncertain by a 
factor of a few. 

Combining the minimum column density of the 
line-of-sight gas with the radial constraints and 
knowledge of the gas having a high covering frac- 
tion yields a mass estimate > 4 x lO^^'Mo for the 
gas emitting the Fe Ka line, assuming the nominal 
radius of emission of 2 x 10^'* cm. The 90 percent 
confidence uncertainty in the distance translates 
into a confidence region on the mass lower limit 
of 5 x 10^^ — 0.1 Mq. If a torus is not the true 
geometry of the gas then, of course, the covering 
factor could be different to this value. Further 
to the uncertainties mentioned, the line may be 
be comprised of contributions from two regions. 
Taking instead the two-component fit to the Fe 
Ka profile then the column and/or global cover- 
ing requirements are reduced for each of the two 
emitting regions. 

The PGA decomposition (jMiller et al. 20091 ) 
suggests a link between the origin of the Fe Ka 
line and that of the line at 5.44 keV. Observation 
of similar lines in other AGN has motivated dis- 
cussion of several possible origins, including spal- 
lation and hotspot emission from the accretion 
disk. 

The solution found in the context of the disk 
hotspot model suggests that the emitting radius 
is between 18 — 24 r^. The orbital timescales at 
these radii are ~ 4 — 6 ks for 18-24 Vg for the black 
hole mass considered here. The observation of a 
steady line flux over 2005 - 2008 provides a con- 
straint on the disk hotspot hypothesis as the three 
year baseline is equivalent to ^ 16, 000 — 24, 000 
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orbits about the black hole over the radial range of 
interest (and tens of orbits just considering the line 
persistence within the 2005 observation). Theo- 
retical modeling indicates that hotspot events are 
not expected to last longer t han a few orbita l 
timescales at these small radii (jKaras et al. 20011 ) 
and that a given hotspot would suffer measurable 
flux and energy chan ges as the material spirals in 
(|Dovciak et al. 20041 ) . However, persistent steady 
lines could arise from a special radius in the disk or 
other rotating reprocessor without any constrain- 
ing expectation of flux or energy variability. If the 
special radius is interpreted as the truncation ra- 
dius of the inner disk this places the innermost 
edge of the disk at 18-24rg. The inner edge of 
the disk may emit more strongly than the rest of 
the structure if it is inflated due to radiation pres- 
sure as in the advection-dominated accretion flow 
scenario. Assuming t he bolometric luminosity to 
be Lboi = lO^^gj-gg-i (jVasudevan fc Fabian 20"09l 
corrected to the TuUy-Fisher distance 15.2 Mpc, 
Russell 20041) and assuming a radiative efficiency 



j^BOL _ Q ^Yie mass accretion rate is estimated 
to be M ~ 0.0035 Moyear-i, ~ 10% of the Ed- 
dington accretion rate. For such high accretion 
rates the transition from thin disk to an advection- 
dominated accretion flo w may occur anywhe re up 
to a radius of ^ 2000rg (jNaravan et al. 19981 ) and 
so our fitted radius is consistent with such a sce- 
nario. If features on the truncated edge of the 
disk are steady in fiux, as found for the line in 
NGC 4051, then Galactic binaries would also be 
expected to show such lines in the low flux state 
and the absence of such would disfavor the trun- 
cated disk origin for the lines. While similar lines 
have not been rep orted to date for Galactic black 
hole binaries (see iDone et al. 20071 for a review) , 
constraints on such lines have not yet been ex- 
plored in that source class, leaving this an open 
question. 

In an alternative model, the fact that the spe- 
cific energy of the line is coincident with Ka emis- 
sion from neutral Cr (5.4 keV) prompts a renewed 
interest in the spallation of Fe as a mechanism for 
enhancing otherwise weak lines, especially since 
PCA is consistent with a common origin for the 
new line and the n eutral component of Fe . In 
a companion paper. Turner fc Miller (2009b( ). ex- 
plore in detail the spallation inter pretation of the 
new result, extending the work of lSkibo (19971 ) in 



the light of new understanding about the environs 
of active nuclei. iTurner fc Miller r2009bl ) find the 
observed abundance enhancement to be high and 
that these extreme enhancement effects are most 
likely to be achieved in gas out of the plane of the 
accretion disk. In such a picture, the timescale 
for spallation may be as short as a few years if the 
cosmic ray output is com parable to the bolometric 
output of the nucleus. iTurner fc Miller (2009bh 
also estimate the expected radio and 7 — ray flux 
from the proposed spallation process in NGC 4051 
and find predictions to be consistent with current 
flux measurements in those bands. 

5. Conclusions 

Our analysis of Suzaku data from NGC 4051 
taken during 2005 and 2008 has revealed line emis- 
sion at 5.44 keV in the rest-frame of the galaxy. 
We have established the reality of the line at 
> 99.9% confidence in data from 2005, supported 
by Monte Carlo simulations that show the proba- 
bility of the line being a statistical fluctuation is 
p < 3.3 X 10~*. The possibility of the line aris- 
ing from a statistical fluctuation in the spectral 
data has been firmly ruled out by establishing its 
detection in time-sliced data. Further to this, we 
have confirmed the line to be evident in all three 
XIS units independently, and established that the 
observed line is inconsistent with arising from the 
X-ray background. 

The source spectrum varies with fiux, and the 
low state is dominated by a hard spectral form 
with the new line plus the neutral component of 
Fe Ka emission superimposed upon that, sugges- 
tive of a common origin for both. The line has an 
equivalent width during 2005 of about 45 eV while 
the Fe Ka line is measured at 195 eV. These re- 
processed signatures show up prominently in the 
source low state when the continuum is suppressed 
by the highest covering fraction of absorption. 

As disk hotspot emission would be expected 
to vary in flux and energy over relatively short 
timescales, the limits on line variability disfavor 
this particular origin although the data remain 
consistent with emission from a special location 
such as the innermost radius of the accretion disk. 
The alternative picture, that the line is Cr i Ka 
emission following spallation of Fe is also found to 
be a good explanation of the data: that possibility 
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and its implications are explored in a companion 
paper. 
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Table 2: The Spectral Model 



Table 1: The Observation Log 



Observation 


Date 


ct s"-^ 


Flux ^ 


Exposure time 


XMM 1 


2001-05-16 


19.27 


2.31 


39.6 


XMM 2 


2002-11-22 


3.15 


0.56 


150.9 


Suzaku 1 


2005-11-10 


0.45 


0.87 


120 


Suzaku 


2005-11-10 


0.32 


0.70 


80.0 


Suzaku l-i2 


2005-11-10 


0.71 


1.26 


39.6 


Suzaku 2 


2008-11-06 


2.06 


2.42 


275 


Suzaku 2-il 


2008-11-06 


1.37 


1.76 


129.2 


Suzaku 2-i2 


2008-11-06 


2.38 


2.75 


92.9 


Suzaku 2-i3 


2008-11-06 


3.63 


3.89 


40.6 


Suzaku 3 


2008-11-23 


1.42 


1.79 


78 



il, i2 etc refer to the intensity selected data, see 
text for details 

^ The observed 2-10 keV flux in units 
10~"^^ergs cm~^ s~^ 

^ Mean exposure time in the CCD instruments in ks 
^ Mean count rate in the 0.5-10 keV band given per 
FI XIS in the case of Suzaku or for the pn in the 
case of XMM 



Parameter 


2005 


2008a 


2008b 


r 


2.47±0.02 


I 


I 


-* ^norm 


0.58 ±0.01 


1.818 ±0.004 


1.283 ±0.005 


"Hi 


2.17 ±0.232 


I 


/ 


Log il 


3.40±0.01 


I 


I 


m2 
" H5 


1.17±0.04 


I 


I 


Log ^5 


n 10-1-0.15 


I 


I 


Covering (Nhs) 


71±2% 


35 ± 2 % 


38 ± 2% 


^ '^■^nbund 


0.40±0.03 


I 


/ 


Per^ 


3.84 ±0.19 


5.41 ±0.26 


4.84 ±0.28 






I 


I 


Log ^6 


2.62±0.05 


I 


I 


Log Reflion^ 


2.94 ±0.07 


I 


I 


Reflionnorm 


1.03 X 10"'^ 


I 


I 



^ PL normalization in units 10 ^photon cm ^s ^ at 1 
keV 

^ Column densities in units of 10^'^ cm^^ 
This zone is ou tflowing at ^ 5500 km/s -see 



Lobban et al. 201Q an d also iCoUinge et al. 2001 



Steenbrugge et al. 20091 



Parameters of the pexrav component, abundance 
as fraction of solar and normalization in units 
lO^^photoncm^^s^^ at 1 keV; 

Normalizations are in photons cm~^ 

I indicates a parameter that was linked to be the same 
for all epochs 

All components are also fully covered by the Galactic 
column of neutral gas with NH{Gal) = 1.3 x lO^'^cm^'^ 
plus three zones of ionized gas determined fr om HETG 
feature measurements (Lobba n et al. 20"lo[ ) and thus 
frozen in this fit. Nhi = 9.6 x lO^^cm^^, log^i = 
0.90; Nh2 = 5.12 x lO^W-^ logC2 = 2.08; Nhs = 
1.08 X 102icni-2,logC3 = 2.14 

Nh5 covers a fraction of the powerlaw (only) 

Nh6 covers the ionized reflector (only) 
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Table 3: X-ray Line Emission in the 3-10 keV band 



Observation 


F6.40 ' 


EW6.40 


'^Xl40 


F5.44 ' 


EW6.44 


AX5.44 


x'/d.o.f. ^ 


XMM 1 




95lil 


143 


0.86^0:86 


23+^'^ 





117/109 


XMM2 


18.07l^:5» 


205lig 


158 


2 50+^-'*^ 
Z.0U_2 09 


07+30 

"^'-30 


3 


146/108 


Suzaku 1 




155_24 


379 


5.0311:1^ 


4611^ 


32 


150/100 


Suzaku 1-il 


IS.SO+i'l? 


201+^^ 


261 


r- 10 + 1. 59 
0-^°-1.60 


46111 


27 


107/100 


Suzaku l-i2 


20.99l;^;;;« 


149li 


127 


4-841I.I? 


26111 


8 


136/100 


Suzaku 2 


no +i.t)b 

Z^.01_j gg 


91+^ 


160 


4 23+ ^ ■'^'^ 


1311 


19 


182/100 


Suzaku 2-il 


21.9512 02 


120+11 


269 


9 or+i.yy 

z.oo_2 00 


lo+l 


4 


135/100 


Suzaku 2-i2 


20.52l|95 


74111 


131 


c; 01 +2.88 

°-°-'--2.88 


16l| 


11 


117/100 


Suzaku 2-i3 


27.751^:11 


75111 


79 


6.8011:1" 


1411S 


5 


114/100 


Suzaku 3 


19.981^:^^ 


107111 


142 


4.311^:^0^5 


17tl^ 


19 


120/100 



il, i2 etc refer to the intensity selected data, see text for details 

^ Line normalization in units 10~^ photons cm~^s~^ for a line having width 

(7 = 50 eV 

^ Fit statistic including both of the lines 

Line equivalent widths are measured against the total observed continuum 
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I I 1 I 1 I 1 I 

4 6 8 10 

Energy (keV) 

Fig. 1. — The offset component from the PCA 
decomposition of all of the 2005 and 2008 Suzaku 
XISO + 3 and PIN data. A strong narrow Fe Ka 
line is evident at 6.4 keV, along with emission at 
5.4 and, more weakly, at 6keV. 




Observed Energy (keV) 

Fig. 2. — Suzaku XISO + 3 data and residuals from 
the mean 2008 Nov 6 spectrum, compared to an 
absorbed powerlaw plus Gaussian line at 6.4 keV 
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Fig. 3. — Suzaku XISO + 2 + 3 data and residuals 
from the mean 2005 spectrum, compared to an 
absorbed powerlaw plus Gaussian line at 6.4 keV, 
showing the residual excess counts at ~ 5.44 and 
5.95 keV 
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Fig. 4. — The very low state data (red crosses) 
from within the 2005 observation shown with the 
model line (solid green), illustrating the emission 
and absorption lines present in the data. 
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Fig. 5.— The 2005 (XIS is black, PIN is red), 2008 
Nov 6 (XIS green, PIN ) and Nov 23 (XIS pale 
blue, PIN magenta) data along with the broad- 
band model described in the text. 
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Fig. 6.— The 2005 data for the three XIS de- 
tectors against the continuum model (black line); 
blue is XISO, red is XIS2 and green is XIS3 
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Fig. 7. — Suzaku data shown as a ratio to the 
mean local continuum model. The red line is the 
data from 2008 Nov 6 and the black line is the 
data from 2005. 
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